Entanglement monogamy and entanglement evolution in multipartite systems 
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We analyze the entanglement distribution and the two-qubit residual entanglement in multipartite systems. 
For a composite system consisting of two cavities interacting with independent reservoirs, it is revealed that the 
entanglement evolution is restricted by an entanglement monogamy relation derived here. Moreover, it is found 
that the initial cavity-cavity entanglement evolves completely to the genuine four-partite cavities-reservoirs 
entanglement in the time interval between the sudden death of cavity-cavity entanglement and the birth of 
reservoir-reservoir entanglement. In addition, we also address the relationship between the genuine block-block 
entanglement form and qubit-block form in the interval. 

PACS numbers: 03.67.Mn, 03.65.Ud, 03.65.Yz 



I. INTRODUCTION 

As an important physical resource, entanglement has 
widely been applied to quantum communication!]] Htl and 
quantum computationiH 0]. It is fundamental to character- 
ize entanglement nature of quantum systems, especially at a 
quantitative level. Until now, although the bipartite entangle- 
ment is well understood in many aspects, the multipartite en- 
tanglement is far from clear ||5|] and thus deserves profound 
understandings. In many-body quantum systems, one of the 
most important properties is that entanglement is monoga- 
mous, which means quantum entanglement can not be freely 
shared among many parties. As quantified by the square of 
the concurrence Jg], a three-qubit monogamy inequality was 
given by Coffman et al. 10] asC^ BC > C\ B +C\ C . Recently, 
its iV-qubit generalization was made by Osborne and Ver- 
straete [8]. Moreover, using some other entanglement mea- 
sures, similar monogamy inequalities have also established 

a Ha eh ee 1 1 411 - However, in these monogamous rela- 
tions, only the single party partition . . . A n is con- 
sidered. Whether it can be generalized to other partitions, such 
as two parties cut AiAj\A^Ai . . . A n , is still an open question 
to be answered. 

On the other hand, the entanglement dynamical behavior 
under the influence of environment is also an important prop- 
erty of quantum systems. This is because that, in realistic 
situations, quantum systems interact unavoidably with the en- 
vironment, and may lose their coherence. It was reported re- 
cently that an entangled state of two qubits interacting respec- 
tively with two local reservoirs would experience disentangle- 
ment in a finite time, even if the coherence is lost asymptot- 
ically lfl5l [l6l ITtI [Hil . This phenomenon is referred to as 
entanglement sudden death (ESD), and has received a lot of 
attentions both theoretically and experimentally (see a review 
paper [ 19] and references therein). 

Recently, Lopez et al. analyzed the entanglement trans- 
fer between two entangled cavity photons and their corre- 
sponding reservoirs, and showed that the entanglement sud- 
den birth (ESB) of reservoir-reservoir subsystem must ha ppe n 
whenever the ESD of cavity-cavity subsystem occurs 12011 . 



However, in this process, whether there exists an entangle- 
ment monogamy relation restricting the dynamical evolution 
is awaited for further studies. Moreover, in the time interval 
where both the cavity-cavity entanglement and the reservoir- 
reservoir entanglement are zero, a subtle issue where the ini- 
tial entanglement really goes is yet to be resolved, although 
the nonzero cavity-reservoir entanglement in this time win- 
dow was pointed out. 

In this paper, based on a new monogamy relation, the entan- 
glement dynamics of two cavities interacting with individual 
reservoirs is studied. It is found that the genuine multipartite 
entanglement is involved in the dynamical process. Particu- 
larly, at a quantitative level, we show the initial cavity-cavity 
entanglement evolves completely to the genuine four-partite 
entanglement in the time interval between the ESD and the 
ESB. In addition, we also address the property of the genuine 
multipartite entanglement which exhibits in the block-block 
form under the bipartite two-qubit partition. 

II. TWO-QUBIT RESIDUAL ENTANGLEMENT AND 
MONOGAMY RELATIONS 

Let us first recapitulate the monogamy inequality in bipar- 
tite single-qubit partition, which can be written as ||8|] 

C 2 A l \A 2 A 3 ---A n - ^AiAi +C\ 1 A 3 + ••■ + < ^A 1 A n - (1) 

The entanglement between subsystems A\ and A2A3 . . . A n 
is quantified by C Al \ A2 A 3 ... An (pAiA 2 A 3 -A n ) 

mmY, x PxTA 1 {PA 1 ), where the r Al (p A J = 2 l l ~ MpIJ 2 ] 
is the linear entropy lull 12211 . and the minimum runs over all 
the pure state decompositions. For the two-qubit quantum 
state pAiAj, its entanglement is analytically expressed as 
C\. A] = [max(0, \/Ai - V^2 - \Aa - VAl)] 2 , with the 
decreasing nonnegative real numbers being the eigenvalues 
of the matrix pij (a y <E> <J y )p*j (o~ y (8 <r y ) [g]. Based on the sum 
of the residual entanglements M Ai = ^a\r(A) ~^2j^A i A j ' 
a multipartite entanglement measure for pure states is 
introduced HHH. 
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Now we analyze the multi-qubit entanglement distribution 
under bipartite two-qubit partition. First, we consider a 2N- 
qubit mixed state pA 1 A' 1 A 2 A' 2 ...A n A' rl with the reduced density 
matrix pAiA'. being a rank-2 quantum state. For this quantum 
state, the following relations hold: 

c 2 

u A 1 A' 1 \A 2 A' 2 ...A n A' 7l 
n 

>^C 2 AiA , ]AiA , (2a) 

i=2 

n n 

^ zJ C ^il^M'» + zL C A' 1 \A i A' i ( 2b > 
i=2 i=2 
n 

> Y.( C MA t +C\ iA , +C\ Ai +C\ A[ ). (2c) 

i=2 

In the derivation of the above inequalities, we have used the 
property that AiA\ is equivalent to a single qubit and the 
monogamy relation in Eq. (1). We here refer to the inequal- 
ities (2a) and (2b) as the strong monogamy relations, and the 
inequality (2c) as the weak monogamy relation. In the rank-2 
case, we define the two-qubit residual entanglement as 

Ma,A' z (PA®«A<®k) = C A,A' t \R(A,A' t ) ~ @) 

where R{AiA'^) denotes the subset of qubits other than 
AiA'i, and i, j in the sum represent the qubit in the sub- 
sets {Ai,^} and {R(AiA' i )}, respectively. It is obvious 
that the residual entanglement is zero when the 2A-qubit 
state is separable under the two-qubit partition. As a non- 
trivial example, we consider the 2A-qubit W state, which 
can be written as \W) 2 n = Oi X \\Q . . .00) + a 2 |01 . . .00) + 
... + a2n|00 . . .01). For this quantum state, we have 

<?MAiWMAn = 4 E-=iEi: 3 K| 2 - K-l 2 and c 2 = 
4|cKi | 2 - \otj\ 2 . Then, according to Eq. (3), the two-qubit resid- 
ual entanglement is zero. Since the square of the concurrence 
is a good entanglement measure for two-qubit quantum state, 
the nonzero residual entanglement M-AiA'. implies the exis- 
tence of multipartite entanglement. 

While for the two-qubit partition of rank-3 and rank-4 
cases, the monogamy relation in Eq.(2) may not hold rf25tl . 

III. ENTANGLEMENT EVOLUTION IN MULTIPARTITE 
CAVITY-RESERVOIR SYSTEMS 

In Ref. lf2(ill . Lopez et al. analyzed the entanglement dy- 
namics of two cavities interacting with independent reser- 
voirs. The initial quantum state of the composite system is 
|$o) = (a|00) +/3|ll)) ClC2 |00) rir2 , where the two entangled 
cavity photons are in a Bell-like state and their correspond- 
ing dissipation reservoirs are in the vacuum states. The inter- 
action Hamiltonian of a single cavity and an A-mode reser- 
voir is H — huja^a + ftJ]feLi ^kb[b k + %Y^k=l 9k{ a b\ + 
bko)). Under the unitary evolution U (H, t) = U Ciri (H, t) ® 
U C2r2 (H, t), the output state is given by 

|$ t ) = a|0000) CiriC2r2 + P\<p t )c iri \<t>t) c 2 r 2 j (4) 



where \<j> t ) = £(t)|10) cr + x(i)|01) cr , and the amplitudes 
= exp(-«£/2) and%(i) = [l-exp(-Ki)] 1 / 2 in the large 
N limit. For this dynamical process, Lopez et al. disclosed 
an intrinsic connection between the ESD of the cavities and 
the ESB of the reservoirs. However, it is not clear whether 
one can establish a quantitative relation of the entanglements 
in different subsystems in the process. Furthermore, it is still 
a subtle issue where the entanglement really goes in the time 
window between the ESD and the ESB. 

We first show that an entanglement monogamy relation 
exists and restricts the dynamical process of the multipar- 
tite systems. The reduced density matrix of a single cavity 
with its reservoir is p Ciri (t) = U Ciri [p Ciri (0)]U^ iri , where 
p cirx (0) = H 2 |00)(00| + |/3| 2 1 10> <10| is a rank-2 two-qubit 
state. Since the unitary operation does not change the rank 
of the matrix, the p Ciri (t) is also a rank-2 density matrix. 
Therefore, the entanglement monogamy relations under the 
partition C\T\\ciTi always hold in the dynamical procedure. 
Particularly, we have 

C 2 Cirilc2r2 (t) >C c 2 iC2 (0+C 2 ir2 W+C c 2 ir2 W+C c 2 2ri (i), (5) 

where the two-qubit entanglements are 

C 2 lC2 (i) = 4[max(|a/?£ 2 |-|/3£x| 2 ,0)] 2 , 
C 2 rir2 (t) = 4[max(|a/? X 2 |-|^ X | 2 ,0)] 2 , 
C 2 Cir2 (t) = C 2 2ri (t) = 4[max(|a/3£xl - |/?£xf , 0)] 2 . (6) 

Here, the bipartite entanglements are quantified by the square 
of the concurrence rather than the concurrence in the analysis 
of Lopez et al. It should be emphasized that, once the ini- 
tial state is given, the bipartite entanglement C 2 iri | C2J , 2 ($t) = 
4|a/3| 2 is invariant in the entanglement evolution, where the 
invariance property of entanglement under local unitary oper- 
ations is used. 

In Ref. rf20ll . the multipartite entanglement is quantified by 
the multipartite concurrence CV c5l . However, Cat is unable 
to characterize completely the genuine multipartite entangle- 
ment. For example, when the quantum state is a tensor product 
of two Bell states, Cat is nonzero. In this paper, we consider 
the two-qubit residual entanglement 

M ciri ($ t ) = C 2 Cirilc2r2 (t) - £C?- (t), (7) 

where i £ {ci,7*i} and j 6 {02,^}. This quantity can not 
only validate the monogamy relation, but also serve as an in- 
dicator of genuine multipartite entanglement in the dynami- 
cal process. According to the expression of |$ t ) in Eq. (4), 
one can deduce that all its three-tangles [7] T^(pijk) = 0, 
because pijk can be written as the mix of a VF-state and a 
product state. Therefore, the nonzero M Ciri (<E> t ) indicates 
only the genuine four-qubit entanglement. In Fig.l, we plot 
the residual entanglement M Ciri as a function of the initial 
state amplitude \a\ and the dissipation time nt. For a given 
value of the a, the M Cin {nt) changes from zero to a maxi- 
mum value, and then decreases asymptotically to zero when 
nit — > 00. Moreover, the maximum values of M Ciri (nt) oc- 
cur in the time nt = ln2 being independent of the amplitude 
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FIG. 1: (Color online) Two-qubit residual entanglement M ciri (4> t ) 
vs the real parameters \ a\ and Kt in the entanglement evolution. 



a. For all possible a, the maximum of the residual entangle- 
ment is M ClTl (a, ln2) = (13\/l3- 19)/34 » 0.81977, where 
|a| = [(9 + VB)/U}^ 2 ~ 0.60889. 

Now, we look into the subtle issue where the initial entan- 
glement goes in the time interval when both cavity-cavity and 
reservoir-reservoir entanglements are zero. We choose the ini- 
tial state parameter a = l/\/T0 which is the same as that 
in Ref. 12011 . and for this value, there is such a time win- 
dow. In Fig. 2, we plot the two-qubit residual entanglement 
M cm and related bipartite concurrences C 2 against the pa- 
rameter Kb. The bipartite entanglement C 2 iri | C2l , 2 (Ki) in the 
process is a conserved quantity (= 0.36) and the monogamy 



relation in Eq. (5) restricts the entanglement evolution. The 
two-qubit residual entanglement M ciri changes from zero to 
the maximum 0.36 in the time [0, — ln(2/3)], then the value 
keeps unchanged until nt = ln3, and finally the M Ciri de- 
creases asymptotically to zero as the time nt — > oo. This 
indicates that the genuine multipartite entanglement is always 
involved in the dynamical process. Particularly, in the plateau 
of Kt € [-ln(2/3),ln3] where all the C?-(t) in Eq. (7) are 
zero, the initial entanglement C 2 iC2 (0) = 0.36 transfers com- 
pletely to the genuine four-qubit entanglement in the com- 
posite system (note that all the three-tangles are zero). In 
this region, the M Ciri is just the C 2 iri \ C2r2 , anc ^ * s ent; angle- 
ment monotone, being able to characterize the genuine four- 
qubit entanglement. For other initial state amplitudes satis- 
fying |a| < \0\/2, there is also a plateau of M Cin («;t) (see 
Fig.l) whose width and value are Ktw = ln(|/3/ct| — 1) and 
M ClTl = A\af3\ 2 , respectively. After a direct comparison, 
we can get that the value is equal to the initial cavity-cavity 
entanglement (C? lC2 (0) = 4|a/3| 2 ) and the width is just the 
time window 112011 between the ESD of cavities and the ESB 
of reservoirs. Here, according to Eq. (6), one can prove fur- 
ther C 2 ir2 (t) = Cl 2Ti (t) = in the interval. Therefore, we 
conclude that the initial entanglement evolves completely to 
the genuine four-partite entanglement in the time window be- 
tween the ESD of cavity subsystem and the ESB of reservoir 
subsystem. 

We also wish to indicate that the nonzero C 2 iri (t) in Fig. 2 
does not come from the initial entanglement C 2 lC2 (0), but is 
generated by a "local" unitary operation U Ciri (H, t) with the 
partition Ciri\c2r 2 . 

IV. BLOCK-BLOCK ENTANGLEMENT VERSUS 
GENUINE MULTIPARTITE ENTANGLEMENT 
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The multi-qubit entanglement property in the plateau region 
is worthy of a further analysis. For the initial state with a = 
1 / a/10, the output state of the evolution can be written as 



i 



|0000) c 



+ 



(8) 



where \ip t ) = ^ (t) j 10) + %(t)|01). Its genuine four-qubit en- 
tanglement is evaluated in bipartite block-block form, i.e., the 
entanglement measure M ciri (V t ) = C\ ir ^ C2r2 = 0.36 
characterizes the genuine block-block entanglement between 
subsystems c\r\ and c^ri- The case for other a with plateau 
region is similar. 

Although the three-tangles Ts(pijk) and the related C 2 in 
the plateau region are zero, the three-qubit subsystems exhibit 
genuine qubit-block entanglements and the relation 



FIG. 2: (Color online) The two-qubit residual entanglement M ciri 
(purple solid line) vs the time evolution parameter Kt, in compari- 
son to bipartite entanglements Cl 1C2 (blue dashed line), C, 2 1? . 2 (green 
dot-dashed line), 10(C;? ir2 + Cl 2ri ) (red dot-dashed line), Cf iri | C2ra 
(black doted line), and Cc iri (yellow solid line) in quantum state 
*t) for which a = l/VWfxJu. 



C 2 

c 1 r 1 \c 2 r 2 



(t)=C 2 Cllc2r2 (t)+C 2 r 



(9) 



holds, in which C c 2 l|c2r2 (i) = 4|a/?| 2 |£(t)| 2 andC 2 i|c2r2 (t) = 

4|a/3| 2 |x(^)| 2 being equivalent to the mixed state one-tangle 
Bg]. This qubit-block entanglement is similar to that of mixed 
states in Refs. 11271 12811 which are entangled but without the 
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FIG. 3: (Color online) The relation between the block-block entan- 
glement and qubit-block entanglement in the plateau region. 

(two-qubit) concurrences and three-tangles. For this kind of 
entanglement, our understanding is that it comes from the gen- 
uine multipartite entanglement in its purified state |28J]. Here, 
Eq. (9) actually presents for the first time a quantitative re- 
lation for understanding the qubit-block entanglement, with a 
schematic diagram being depicted in Fig. 3. 

V. DISCUSSION AND CONCLUSION 

Entanglement monogamy is a fundamental property of mul- 
tipartite entangled states. We argue that the violation of the 
monogamy relations in Eq. (2) for higher rank cases is be- 
cause the square of the concurrence does not have the additiv- 
ity, i.e., C 2 AiA , ilMA , 2 ? C\ iM + C\, iA , 2 for the tensor prod- 
uct of two Bell states. The von Neumann entropy has this 
additivity property, however, it has the negative residual en- 
tanglement for multipartite systems ll29ll . How to define an 
additive entanglement measure with nonnegative residual en- 
tanglement is still challenging. 

The monogamy relations in Eq. (2) can be applied to other 
systems lf30h only if the individual system-environment is in a 



rank-2 quantum state and the evolution has a tensor structure 
U(H, t) = U SlEl {H, t) ® U S2E2 (H, t) ® . . . ® U SnEn (H, t). 
Moreover, based on this relation, one can derive other useful 
monogamy inequality. For example, if the initial state of a 
three cavity-reservoir composite system is \^q) = (a|000) + 
/3|lll)) c |000) r and the individual cavity-reservoir interaction 
is the same as the previous one, we can derive 

C c 1 r 1 | C2 r 2 c 3 r 3 ( ) ^ T 3 (Pc lC2 c 3 (*) ) + T 3 (p ri r 2 r 3 (*) ) , (10) 

where ^^^^(O) = 4|a/3| 2 gives an upper bound for 
the three-tangles in the entanglement evolution. 

In conclusion, we show that a monogamy relation restricts 
the entanglement evolution of two cavities with individual 
reservoirs. Moreover, based on the relation, we find the initial 
state entanglement evolves completely to the genuine four- 
partite entanglement in the time interval between the ESD 
of cavity-cavity entanglement and the ESB of the reservoir- 
reservoir entanglement. In addition, we give a quantitative 
relation between the block-block entanglement and the qubit- 
block entanglement. 
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